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Inclusive electromagnetic reactions in few— nucleon systems are studied basing on 
accurate three- and four— body calculations. The longitudinal *He(e,e') response 
function obtained at g < 600 MeV/c completely agrees with experiment. The 
exact *He spectral function obtained in a semirealistic potential model is presented, 
and the accuracy of the quasielastic response calculated with its help is assessed, 
as well as the accuracy of some simpler approximations for the response. The 
photodisintegration cross section of ^He obtained with the realistic AV14 NN force 
plus UrbanaVIII NNN force agrees with experiment. It is shown that this cross 
section is very sensitive to underlying nuclear dynamics in the ~ 100 MeV 
region. In particular, the NNN nuclear force clearly manifests itself in this region. 



1 Introduction 

Our recent results on the '*He(e, e') process and the photodisintegration cross 
sections of trinucleons are presented in the talk. The issues of interest are 
the accuracy of the conventional 4N formulation for the ''He(e, e') reaction 
and the accuracy of its description in the framework of the spectral function 
approximation, the sensitivity of the photodisintegration cross section to the 
choice of NN interaction, and the manifestation of NNN force in the process. 
The quantities we study are response functions having e.g. in the (e, e') case 
the form 

Riq,u;) = J df\{^f\diq)\^oiMo))\'6[u; - (Ef + / {2Am) - E^)], (1) 

q and oj being momentum and energy transferred to the nucleus. We avoid 
calculating the complete set of final states ^ff^j^nd we obtain the response 
functions with the method explained in refs. ufl requiring only bound-state 
type calculations. Within the adopted formulation of the problem the numeri- 
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cal results below are accurate at the per cent level. 



2 The (e, e') response and spectral function of '*He 

In the present report, we present the results on the '*He(e,e') respoiise. in. the 
q range of 300-600 MeV/c. Additional information can be found in BB'QO. A 
non-relativistic 4N Hamiltonian with central NN potentials reproducing the 
^Sq and ^5*1 NN phase shifts up to high energy hpii been used as a dynamic 
input. This Trento (TN) potential is listed in Ref □. The conventional single- 
particle transition operator that includes the Sachs form factors G^j" with 
the usual relativistic correction has been adopted. No further approximations 
have been used, and the FSI is thus fully taken into account. In Fig. 1 the 
calculated longitudinal response function cj) is shown with solid curves 

r four values of the transferred momentum q [MeV/c] along with the Bates 
and Sajday Q data. (The results for q < 500 MeV/c differ^ bit from those in 
Refs. due to minor corrections in the formulation, see □.) One observes a 
very good agreement with experiment in the whole range of q values. (A drop 
in the data at g = 500 MeV/c in the vicinity of lo =250 MeV is at variance 
with the sum rule systematics.) 

To interpret the results obtained let us estimate the maximal momenta 
in the initial state contributing substantially to the response. Let us con- 
sider e.g. the final state component in which only the N-(yl — 1) relative 
motion is fast. It is sufficient here to consider this relative motion in the 
plane wave approximation. Let k^ez be the relative momentum between the 
fast nucleon and the remaining — 1) subsystem. The FSI admixes the 
momenta lying in the range k^ez + Ak to k^ei. The momenta Ak are de- 
termined by the nuclear force, and they are lower than ^^rel values of interest 
for us. We have Kei = [{A - l)M][kAr - {A - 1)-'^J^a-i] k^ - ^-^K^. 
The corresponding single-particle momentum probed in the ground state is 
kj„iss = krei — — ^) I A. The predominant contribution comes from the 
states with k^e; directed along q. Using the energy conservation in the form 
u) ~ kl^iA/[2[A - l)m] + q^l{2Am) - Eg^ and setting uj + Eg^^ q'^/{2m) + A 
one obtains k^i^s ^ [{A - 1)/A]\{q'^ + [A/{A - l)]2mA}i/2 - q\. Let us con- 
sider the right wing of the response. Setting, for example, q = 600 MeV/c, 
A =100 MeV, c.f Fig. 1, one gets kmiss — 0.7 fm~^. As a consequence, only 
low momenta in the nuclear interaction play a role in the kinematical region 
considered. Contributions from other reaction mechanisms cannot change this 
conclusion. The semirealistic NN potentials of the type we used probably lead 
to the same low momentum content of the wave functions as the fully realistic 
NN interactions. This is illustrated in Fig. 2 where the single-particle momen- 
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Figure 1. Longitudinal response of ^He in MeV~^. Full results: solid curves, PWIA calcula- 
tion: dashed curves. Bates and Saclay data are represented with open diamonds and closed 
circles, respectively. 



Figure 2. Total and partial momentum distributions of ^He with the TN potential and with 
Argonne uig+Urbana IX . 



turn distribution of *He for opj NN potential is compared with that obtained 
from a rcahstic nuclear force EJ. The agreement above with experiment at the 
right wing of the response then testifies to a correct description of the corre- 
sponding low-momentum content of the wave functions within the adopted 4N 
formulation. It shows the applicability of the conventional form of the electro- 
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magnetic interaction as well. On the contrary, rather high nucleon momenta 
in the ground-state wa.ve function are probed at the left wing of the response 
at high q and small u; a. 

The SF of '''He and the corresponding QE response have also been cal- 
culated exactly with the same NN potential. The SF S{k,E) represents the 
joint probability of finding a nucleon with momentum k and a residual {A — 1) 
system with energy E + Egs : 

StAk,E)^Y.Mo,sJ df\{^f-';ks,^^^iMo))\'S[E-{Ef-'-E^)]. (2) 

Here and are the spin and isospin quantum numbers of the nucleon, E^~^ 
and '^f~^ are eigenvalues and eigenstates of the {A ~ 1) system, and E^ and 
are such quantities for the ground state. The energy E can be viewed as 
the removal energy of the nucleon. The calculation of the SF from its definition 
requires obtaining the complete set of states of the residual (A—l) subsystem. 
Previously the spectral function has been calculated in the A=3 case, which 
is essentially the two-body problem, and, approximately, for nuclear matter, 
see □ for the references. We have avoided calculating the complete set of 
A — 3 continuum states with the help of the approach mentioned in Sec. 1. 
The proton SF Sp{k,E) thus obtained is shown in Fig. 3. It includes two 




Figure 3. SF Sp{k,E) of *He with TN potential in units of fm^ McV"^ 

contributions: from the A — 3 rest subsystem in the ground state, and in 
the continuum, c.f. Eq. below. The first of these is determined by the 
corresponding momentum distribution, and the second one is shown in the 
figure. The values of Sp{k,E) > E^^^ + 1 MeV are plotted. (We note that 
S'p(fc, E^^) = and thus Sp{k, E) exhibits a strong slope at low energy.) For 
momenta below 2 fm^^ one finds a sharp maximum at about 2 MeV above 
Efj~^ . In Fig. 4 the cut of the calculated spectral function at fc = 0.25 
fm~^ is shown. These results appear to have subsequently been confirmed 
experimentally Ell. For k > 2 fm^^ the calculated spectral function exhibits a 
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Figure 4. The cut of the "He SF at fc = 0.25 im'^. The values of SF are in IQ-^- fm^ 
MeV-i. 



ridge where the peak position shifts to higher E for increasing k. This ridge 
should presumably correspond-to a proton struck from the region of a strong 
two-nucleon correlation, see BO. 
The QE response 



is shown in Fig. 1 with the dashed curve. At q 600 MeV/c the height of the 
peak is overestimated by about 9% in the framework of the QE approximation. 
At low uj, the PWIA assumption inherent to R'j^ obviously breaks down. For 
high UJ, the relative role of disregarded lower momenta admixtures in the final 
state wave function increases with to. Due to these reasons the relative error 
of the QE approximation (^) increases at the wings of the response. It is seen 
from Fig. 1 that the QE peak is shifted to higher energies with respect to the 
true one. The shift of the peak can be qualitatively understood considering a 
nucleon at rest in a potential well: uj can be estimated as g^/ (2m) + Vf{q) — Vi, 
where Vf^i are the potential energies before or after interaction with the virtual 
photon. While Vf is negative, it becomes zero in PWIA leading to an increase 
in UJ. 

It is advantageous to have a simple and good approximation for the QE 
response. One can write 



Here the first term represents the residual A — 3 system in the ground state, 
and we do not differentiate between the t-p and n-'^He channels. The second 
term is depicted in Fig. 3 and represents the contribution of the residual A = 3 
system being in continuum. It is seen from Fig. 3 that for low k relevant 




(3) 



S{k, E) = ntp{k)5[E - E^i^K) + £;o(^He)] + S,nei{k, E). 



(4) 
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to us almost all the strength in Smei is concentrated close to the thershold 
E =^ Emin- This suggests the approximation (to be put into Eq. (^) 

S{k, E) ^ nfp{k)5[E - ^o('H) + E^^c)] + nt>j,{k)5[E - Et>p + Eo^lle)] (5) 

where Efp is the threshold breakup energy of the residual nucleus. Using (||) 
one does not need continuum state calculations since nfp{k) = n{k) — ntp{k). 
One obtains an even simpler approximation considering that Et*p ~ -Bo("^H): 

Sik, E) ~ n{k)5[E - So(3h) + ^oC^He)]. (6) 

Eq. (0) was used in the literature (see e.g. 0). One may note that at the small 
k values of interest the term Sinei m (^ is only of secondary importance (e.g. 
the integrated values of ntp and rifp equal 0.89 and 0.11, respectively,) which 
improves the quality of the approximations (j^), (^). In Fig. 5 we show the 




Figure 5. QE response (^) with the SF of Eq. (|5|) (dashed curve) and of Eq. (|6|) (dotted 
curve) relative to that with the full SF (the QE peak is marked by an arrow). 



responses (g) obtained with these approximations relative to Eq. (^) with the 
full spectral function. Except for small lo where the wrong threshold behavior 
in Eq. (^ may play a role, the three responses are very similar, particularly 
in the peak region. (At q = 300, 400, and 1000 MeV/c one has very similar 
results.) 

3 Total photodisintegration cross sections of trinucleons 



Earlier theoretical work on the photodisintegration cross section of trinucleons 
in a wide energy range is outlined in Ref. L3. The first rigorous calculations 
have been performed with separable central potentials. Mixed symmetry states 
have been disregarded even though they play an important role in the process. 
Approximate calculations of the three-body breakup with low-energy central 
local potentials were also performed. In Ej the cross section has been obtained 
with a non-central potential but only the three lowest hyperspherical terms in 
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the expansion oLcontinuum wave functions were retained which may influence 
the results. In Il3 the '^H two-body photodisintegration differential cross sec- 
tions at 9 = 90° wei£ calculated with Paris and Bonn NN realistic potentials for 
Ey < 40 MeV. In llj accurate calculations of the totaLcross sections have been 
performed with central TN and Malfliet-Tjon (MT) til potentials reproducing 
the s-wave phase shifts up to high energy 

In the present work we calculate the total photodisintegration cross section 
with realistic NN irtteractions for the first time and we include an NNN force 
in addition. As in til, the approach outlined in Sec. 1 has been used. The 
calculations are performed for photon energies up to pion threshold, and we 
pay special attention to the high energy tail of the cross section. We obtain a 
sizeable contribution of NNN force to the cross section and also comment on 
the sensitivity to the NN force. It may be natural to expect some sensitivity 
to the NN force since if one takes, for example, =110 MeV and considers 
the N+d breakup then the N—d relative momentum is about 1.9 fm~^, and 
single-particle momenta in 0^*0 of this size may be probed in the process. 

The calculation was done with a non-relativistic three^ucleon Hamilto- 
nian. The realistic NN interactions AV14 E£l and TRSB lla were employed. 
These NN interactions underbind "^H by about 0.8 MeV and 1 MeV, respec- 
tively. In the third version of the calculation, the AV14 NN interaction plus. 
NNN UrbanaVIII force are used. Along with these results the cross sections til 
with central TN and MT potentials are listed as well. 

For the latter potentials, as well as in the case of the softer-core TRSB 
interaction, we used the Jastrow correlated hyperspherical basis c2l. The cor- 
relations speed up the convergence but not to a sufficient degree in case of the 
stronger core AV14 interaction. In the Jastrow case, one chooses the two-body 
correlation function as the mean value of the triplet and singlet NN low-energy 
wave functions. However, for interactions like AV14 it is advantageous to take 
into account differences between correlations in different two-nucleon spin- 
isospin states. One possibility is to use different Jastrow factors at different 
groups of basis states before symmetrization Eil. However, after symmetriza- 
tion the correlations then lose their role as correlations in given NN states and 
become some_fifFective correlations which still may substantially speed up the 
convergence Ed due to the many variational parameters involved. Such a pro- 
cedure is not convenient in our case since our response calculations are not 
variational ones. In the present work, in the AV14 case, we use a new version 
of a correlated basis. Spin-isospin dependent correlations are used, and they 
retain their role as correlations in the given four spin-isospin NN states, see 
LJ. In Table 1 the properties of our ground state wave function of '^H are listed 
in comparison with Ref. CJ. The UrbanaVIII potential was taken as the NNN 
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Table 1. Bound-state properties of '^H with the AV14 interaction 





Et, MeV 


frinsf fill 


P(D), % 


£;i,{NN+NNN), MeV 


present rwprk 


7.71 


1.66 


8.85 


8.52 


Ref. y 


7.67 


1.67 


8.96 


8.46 



force in the last column. 

In the present work we confine ourselves with the unretarded El approxi- 
mation for the transition operator which is known to be a very good one for the 
total cross section at least for not too high energies. Within this approximation 
the transition current operator may be rewritten exactly as the dipole moment 
operator of a system, while MEC currents are taken into account automati- 
cally. In Fig. 6 our total photodisintegration cross section of ^He calculated 



■■■■f i- 







Figure 6. Photodisintegration cross section of ^He. Theoretical curve is obtained with the 
AV14+UrbanaVIII nuclear force (full curve). Experimental data are of Refs. (filled circles) 
and (upper and lower bound delineated with dotted curves). 



with AV14+UrbanaVIII interaction is shown along with experimental data. 
The theory leads to an agreement with the data of Ref. Eil. In Fig. 7 our 
results in the tail region of the "^H cross section are shown. Separate cross 
sections for transitions to final states with isospin T = 1/2 and 3/2 are pre- 
sented. Transitions to T = 3/2 correspond to the three-body breakup while 
those to T = 1/2 are known to correspond predominantly to the two-body 
breakup. One sees a large difference between predictions of different models 
in this region. In the T = 1/2 case, realistic NN interactions lead to much 
higher cross section than central NN forces. This can be attributed to tensor 
correlations produced by realistic interactions. It is also seen that the NNN 
nuclear force clearly manifests itself by leading to a sizable further increase in 
the cross section. As to the T = 3/2 case, the interesting point is that the cross 
section obtained with the realistic TRSB NN interaction is close to those for 
the central NN potential models but strongly differs from that obtained with 



tri99: submitted to World Scientific on February 9, 2008 8 



0.1 o 
0.08 



0.04 - 



AV14+NNN 

AVI 4 

TRSB 

TN 

MT 



AV1 4+NNN 

AV1 4 

TRSB 

TN 

MT 



80 100 1 20 SO 100 120 140 

to LJVleVJ CO {IVleVJ 

Figure 7. Photodisintcgration of •^H. Results with various choices of nuclear interactions. 
The UrbanaVIII potential was taken as the NNN force. 



AV14. We are planning to consider other realistic NN interactions as well to 
clear out a sensitivity to the choice of NN interaction. 

There exist only two experimental data points in the energy range pre- 
sented in Fig. 7 (see Fig. 6, the ^He and cross sections are practically 
indistinguishable in this region EJ) . These experimental data are in favour of 
the AV14-|-UrbanaVIII interaction. Our results thus show that, in contrast 
to many quantities conventionally studied in few-nucleon physics, an experi- 
mental study of the process in the considered energy region may allow a clear 
discrimination between various models of NN interaction and detection of the 
contribution of NNN force. Corrections from retardation and the contribu- 
tion of other multipoles will be studied separately. They are rather small, and 
they cannot change the conclusion about high sensitivity of the process to nu- 
clear dynamics. Obviously, one can anticipate even more sensitivity to nuclear 
dynamics in exclusive experiment. 
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